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We review recent developments for the description of far-from-equilibrium dynam- 
ics of quantum fields and subsequent thermalization. 



1. Far-from-equilibrium quantum fields 

The advent of relativistic heavy-ion colhsion experiments as well as rapid 
progress in cosmological observations triggered an enormous increase of 
interest in the dynamics of quantum fields out of equilibrium. Much 
progress has been achieved for systems close to thermal equilibrium us- 
ing (non-)linear response techniques or with effective descriptions based on 
a separation of scales in the weak coupling limit^'^. These methods provide 
an efficient description of the dynamics on sufficiently long time scales in 
their range of applicability. Major open questions for our theoretical un- 
derstanding are currently highlighted by experimental indications of early 
thermalization in collision experiments starting from extreme nonequilib- 
rium situations'^'^. They provide an important motivation to study the 
far-from-equilibrium dynamics of quantum fields and subsequent thermal- 
ization from first principles. 

There are few necessary ingredients to describe far-from-equilibrium 
quantum fields. Firstly, in contrast to close-to-equilibrium field theory, the 
initial density matrix can deviate substantially from thermal equilibrium. 
Fully equivalent to the specification of a density matrix is a description 
of the noncquilibrium initial conditions in terms of correlation functions. 
Secondly, noncquilibrium time evolution involves no other dynamics than 
the one dictated by the underlying quantum field theory. In particular, the 
dynamics can be conveniently obtained from the effective action F for given 
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initial correlation functions. Since F is a Lcgendre transform of the gen- 
erating functional for connected Green's functions, there are very efficient 
approximation schemes available. 

Practicable and systematic approximations for nonequilibrium situa- 
tions may be based on the two-particle irreducible (2PI) effective action'"''^'^. 
It has been recently demonstrated for scalar^'^'^" and fermionic^^ quantum 
field theories that a systematic coupling-expansion^ or l/A''-expansion^'®'^^ 
of the 2PI effective action can describe far-from-equilibrium dynamics and 
subsequent thermalization without further assumptions. At lowest nontriv- 
ial order these approximations contain scattering contributions including 
off-shell and memory effects. In particular, the 2PI l/A''-expansion has 
been shown to solve the problem of an analytic description of the dynamics 
at nonperturbatively large densities^^'^^. 

A related approach that has been successfully applied^^ to nonequilib- 
rium quantum field theory is motivated by truncating Schwinger-Dyson 
equations^^. It is equivalent to an approximation including the complete 
next-to-leading order (NLO) contribution and part of the NNLO graphs 
of the 2PI l/A'^-expansion^^. Similar to the 2PI effective action technique, 
promising approximation schemes may also be obtained in a systematic way 
by considering so-called two-particle "point-irreducible" (2PPI) graphs. In 
the context of nonequilibrium this has been investigated in Ref. 16. 

These methods provide a systematic and practical tool to go beyond 
standard mean-field type approximations^^, such as Ilartree(-Fock) or lead- 
ing order large-iV. Although the latter present a useful guide for various 
problems, they e.g. fail to describe exponential damping of correlations at 
early times or the late-time dynamics leading to thermalization^. These 
problems are related to the appearance of an infinite number of conserved 
quantities, which are not present in the fully interacting or finite-Ar field 
theory. The spurious constants of motion restrict the nonequilibrium evo- 
lution at all times and prevent an approach to thermal equilibrium. Al- 
though a straightforward task in principle, going beyond mean-field type 
approximations has long been a major difficulty in practice: Similar to per- 
turbation theory, standard approximations based on 1/A'^-expansions of the 
one-particle irreducible (IPI) effective action can be secular in time. Even 
for small couplings these approximations break down after a short time and 
do not describe thermalization^^. In the following we review two-particle 
irreducible expansion schemes and recent nonequilibrium applications. 



*For somewhat improved results using inhomogeneous mean fields see Ref. 18. 
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2. 2PI effective action for scalar and fermionic fields 

Consider a theory describing N scalar fields cpa and Nf fermion flavours V'i 
with classical action {a,b= 1 . . .N; i,j = 1 . . .Nf) 

S[ip;i;,i;]= J d^x[^d^^ad''<Pa + Mx)i^Mx) + V{<p;iP,i;)} . (1) 

Here V{ip; -ip, ijj) parametrizes the interaction part to be specified later, and 
^ = 7^9'^ with Dirac matrices 7^ (/i = . . . 3). All Green's functions of this 
theory can be obtained from the 2PI effective action^, which is parametrized 
by the connected one- and two-point functions: 

(j>a{x) = {ifia(x)) , Gab(x,y) = {T^Lpa{x)!fib(y)) - (l)a(x)(f>b(y) , 

and similarly for the fermionic fields. Here, the field operators are time- 
ordered along a closed time path C, as depicted in Fig. 1. For a vanishing 
fermionic "background" field the 2PI effective action can be written as^ 

r[.^, G; D] = S[^; 0, 0] + ^Tr, In G'^ + ^Tt, Go 

-iTr^lnD-^ -iTi^DQ^D + T2[(l>,G;D], (3) 
where the classical inverse propagator for the scalars is given by 



(4) 

and equivalently for the classical inverse fermion propagator iD^^. The 
term T2[(t>, G; D] in Eq. (3) contains all contributions beyond one-loop or- 
der and can be represented as a sum over closed 2PI graphs only^, i.e. dia- 
grams which cannot be disconnected by opening two propagator lines. The 
equations of motion are obtained by extremizing the effective action: 

6mG;D] 5mG:D] (5r[0,G;D] 

S(f)a{x) ' SGab{x,y) ' SDij{x,y) 

The resulting equations of motion can be written as differential equations 
suitable for initial value problems (for details see e.g. Ref. 9). All higher 
n-point functions can then be obtained from the generating functional 
T[(I),G;D] at any time during the evolution. Without further truncation 
the effective action contains the complete information. 
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Figure 1. Traces involve a finite closed 
time path C suitable for the initial value 
problem with a causal time evolution. 
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3. Systematic approximation schemes 

A major advantage of the 2PI effective action is that it aUows one to ob- 
tain suitable approximation schemes for nonequihbrium problems in a sys- 
tematic way. As an example, we consider first a scalar field theory with 
0(A'')-symmetric classical action given by Eq. (1) with 

yif) = ^ fa^a + {fa<Paf ■ (6) 

One can classify the various contributions to the 2PI effective action ac- 
cording to their scaling with N . More generally, one can write the sum of 
2PI diagrams as 

r2[0, G] = tY'[4>, g] + r^LO[</), G] + r^LOf^:^ ^ ^ ^ ^ 

where T^o ~ N, r^^o _ ^nnlo _ ^^jy g^^j^ ^^^^^ ^^^i 

following contribution is suppressed by an additional power of 1/A''. The 
leading order contribution corresponds to only one diagram and reads^ 

^2''[G]=-^ J^G,a{x,x)G,,{x,x). (8) 

Here we use the notation = J d.x" J dx. The full ncxt-to-lcading contri- 
bution contains an infinite series of 2PI diagrams which can be summed:*'^^ 




Figure 2. The complete NLO contribution to r2. The dots indicate that each diagram 
is obtained from the previous one by adding another "rung" with two propagator lines 
at each vertex. The crosses denote field insertions. 
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and the function l{x,y;G) sums an infinite series of "rungs" (cf. Fig. 2): 

A X f 

l{x,y;G) = —Gab{x,y)Gab{x,y) - i— J I{x,z;G)Gab{z,y)Gab{z,y) . 

The 2PI l/7V-expansion provides a controlled nonperturbative approxima- 
tion scheme. It can be in particular applied to critical phenomena near 
continuous phase transitions, or to nonpcrturbativcly large densities as de- 
scribed below^^'^^. Instead, for sufficiently dilute systems and weak cou- 
plings a 2PI coupling-expansion can be appropriate. The lowest non-trivial 
order then contains all diagrams with topology as given by the first graph 
in the upper and lower series presented in Fig. 2. 

The systematic expansion schemes discussed above can be straightfor- 
wardly applied to fermionic theories as well^'^^. As an example, consider a 
theory with two fcrmion flavors coupled to an 0(4)-vcctor of scalar fields 
ipa = (c, 7?) via a chirally invariant interaction with Yukawa coupling g: 

V"(v5;V7,-0) = (o-^ +7r^) ~ g'tp[(7 + ij^f ■ T:]ip . (12) 

This is a simplified version of the linear sigma model, with no scalar self- 
interaction, which is sufficient to obtain thermalization as is described 
below^^. The 2PI effective action at lowest nontrivial order in a systematic 
coupling expansion then contains only the two-loop graph shown in Fig. 3. 



Figures. Two loop contribution to 

I .] The solid and dashed lines represent the 

V J full fermionic (D) and bosonic (G) prop- 

agators, respectively. 

For discussing nonequilibrium results it is useful to rewrite the propagators 
in terms of the spectral function p and the statistical function F'^^'^: 

G{x, y) = F{x, y) - '-p{x, y) sign^ {x° - y") (13) 

and equivalently for the fermion propagator D. We emphasize that out of 
equilibrium F and p are independent functions. However, in thermal equi- 
librium both functions are related by the fluctuation-dissipation relation, 
which can be written in Fourier space^°'^: 

F(-^)(a;,p) = -i(n^^{w) + p(^^)(a;,p-) , (14) 

where nBE(w) = ie'^^'^ — denotes the Bosc-Einstein distribution with 
temperature T. Starting far from equilibrium, Eq. (14) can be used to test 
thermalization at sufficiently late times^^ as is discussed below. 
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4. Comparison of LO, NLO cind exact results 

LO vs. NLO: Why is it crucial to go beyond leading-order large- TV or 
Hartree approximations? We consider for a moment the LO contribution, 
Eq. (8), or the Hartree approximation which also includes the first graph 
of Fig. 2. This takes into account space-time dependent mass corrections 
and neglects, in particular, direct scattering. An important consequence is 
the appearance of an infinite number of conserved quantities, which are not 
present in the fully interacting or finite- AT theory. These can be written 
as9'i7 

F{t,t';p)dtdt'F{t,t';p) - {dtF{t,t';p)fy^^ = no{p) + J , (15) 

-I t=t' 2 

where we have considered a spatially homogeneous situation with F{t, t';p) 
denoting the spatial Fourier transform. These additional constants of mo- 
tion can have a substantial impact on the time evolution, since they strongly 
constrain the allowed dynamics. As a characteristic example, one may con- 
sider the time evolution for so-called "tsunami" initial conditions with a 
vanishing field expectation value: A large initial particle number density 
in a narrow range around a characteristic momentum p = ±pts- Here we 
consider first 1-1-1 dimensions, where the initial condition is reminiscent 
of two colliding "wave packets" with opposite and equal momentum. The 
left graph of Fig. 4 shows the equal-time two-point fmiction F{t,t:p) as a 
function of time both for the LO and NLO approximation^. One observes 
that at NLO the highly populated modes around pts get quickly depopu- 
lated, while the density for low momentum modes increases. In contrast. 




Figure 4. Comparison of results at LO and at NLO in the 2PI l/AT-expansion (iV = 10). 
Left: Equal-time two-point function F{t, t;p) for different momenta as a function of time 
for "tsunami" initial conditions. Right: Unequal-time two-point iunction F{t,0;p = 0) 
as a function of time following a sudden "quench" . (All in units of the initial-time mass.) 
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at LO the evolution of the equal-time correlator is trivial. This reflects the 
fact that at LO or for Hartree approximations, there exists a well-defined 
particle number for each momentum mode, which does not change in time:^ 
no{p) of Eq. (15). This is in sharp contrast to the NLO approximation or 
the fully interacting real scalar theory, where there is no conserved particle 
number and the LHS of (15) depends on time. Similarly strong qualitative 
differences can also be observed for unequal-time correlation functions. As 
an example, we show on the right of Fig. 4 the two-point correlation with 
the initial time, F{t,0;p = 0), following a sudden "quench"^. After some 
oscillations one finds that at NLO the unequal-time correlator quickly ap- 
proaches an exponentially damped behavior. In contrast, at LO no effective 
suppression is observed. This reflects the fact that at LO the presence of 
additional constants of motion keep the details about the initial conditions. 

NLO vs. exact results: In view of the substantial changes encoun- 
tered by going from LO to NLO for finite A'', one wonders what happens 
at NNLO or beyond. One can rigorously answer this question in a well- 
defined limit: Nonequilibrium dynamics in classical statistical field theory 
can be solved exactly up to controlled statistical errors, using numerical in- 
tegration and Monte Carlo techniques^^'^*'^^'^^'^^. The 2PI l/A'-expansion 
can be equally well implemented in the classical as in the quantum casc^^. 
Therefore, in the classical field limit one can compare with results includ- 
ing all orders in 1/N. In particular, for increasing occupation numbers per 
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Comparison between exact (MC) results and results from the 2PI 1/N— 
expansion at NLO both for the classical and the quantum field theory. Left: Time 
evolution of the equal-time correlation function for an initial condition similaj: as for 
Fig. 4 with high occupation numbers {N = 4). Right: Damping rate, extracted from 
the unequal-time correlation function (cf. Fig. 4, right). Here the initial conditions are 
characterized by low occupation numbers so that quantum effects become sizeable. 
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mode the classical and the quantum evolution can be shown to approach 
each other, if the same initial conditions are applied^'*: For sufficiently high 
particle number densities one can strictly verify how rapidly the 1/N series 
converges. 

Fig. 4 shows a comparison of exact and NLO rcsults^^. One observes 
that the NLO result follows closely the exact solution already for moderate 
values of N. This concerns both the far-from-equilibrium dynamics at early 
times as well as the late-time behavior close to equilibrium. The 2PI 1 /N- 
expansion exhibits very good accuracy already in lowest nontrivial order! 
We emphasize that also for low densities, where the classical field approxi- 
mation provides not a good description for the quantum theory (cf. Fig. 4), 
a coupling-expansion to three- loop and a 1/iV-expansion to NLO give very 
similar results for sufficiently weak coupling and not too small N >2. One 
is lead to conclude that, once the spurious conserved quantities of LO or 
Hartree approximations arc removed by interactions, the results are rather 
independent of the details of the approximation. 

1/N -expansion in the limit N = 1: The agreement between NLO 
and exact results is found to be rather good for values as small as > 2 
(cf. Fig. 5)^^. However, in the limit N = 1 subleading contributions in 1/A'' 

are no longer suppressed and, in particular, the 1/A'^-cxpansion at NLO ex- 
hibits a reduced damping rate*^. The case N = 1 has been investigated in 
classical field theory in Ref. 25 for initial conditions with zero and nonzero 
field (p. An approximation, which includes part of the NNLO corrections, 
such as BVA^^, is found to give a better damping rate for the field. How- 
ever, a worse result is obtained for the two-point function, which reflects the 
absence of an expansion parameter for AT = 1. For weak interactions a sys- 
tematic expansion can be based on the 2PI coupling-expansion for N = 1. 

Spontaneous symmetry breaking: It is a well-known fact that there 
is no spontaneous symmetry breaking in one spatial dimension at nonzero 
temperature or energy density. In contrast, LO or Hartree approximations 
can exhibit a phase transition independent of the number of dimensions. For 
nonequilibrium. Cooper at al. have pointed out in Ref. 10 that the 2PI 
expansion at NLO shows no spontaneous symmetry breaking in one spatial 



''We note that this can be understood from the fact that the three-loop contribution to 
the 2PI effective action (which gives the lowest order contribution to the damping rate 
in a coupling expansion) is ~ {X/N)'^[l + 2/N]. Comparing this with the corresponding 
three-loop contribution at NLO in the 1/iV-expansion^'^^ for N = 1, one finds a factor 
of three less for the latter (cf. Sect. 3). 
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Figure 6. Examples of nonequilibrium time evolutions of the field expectation value </>(t) 
for one (left^") and for three (right) spatial dimensions using different approximations. 
The 2PI l/AT-expansion (here for the limit N = 1) correctly describes the absence of a 
spontaneously broken phase in one spatial dimension, with ip{t) — > at late times. In 
contrast, spontaneous symmetry breaJfing can be observed at NLO in 3 -f 1 dimensions. 



dimension, thus curing the problem of mean-field type approximations. As 
an example, the left graph of Fig. 6 (taken from Rcf. 10) shows their results 
for the time evolution of the field expectation value (pit). One observes 
that, while the Hartree approximation oscillates around a nonzero field 
value, the NLO result tends to the symmetric phase with = 0. We 
emphasize, however, that the latter results have been obtained from a 1/A'^- 
expansion in the uncontrolled limit A'^ = 1. This becomes obvious once e.g. 
a part but not all of the NNLO contributions are taken into account, as 
denoted by BVA in Fig. 6, where again one finds spontaneous symmetry 
breaking in one spatial dimension. More recent results, demonstrating the 
absence of a nonequilibrium phase transition in one spatial dimension have 
been obtained from a loop-expansion in a related (2PPI) approximation 
scheme^^. It is interesting that the systematic 1/A-expansion is found to 
give qualitatively correct results even for Af = 1. Here we show for the first 
time that the 2PI 1/A-expansion exhibits spontaneous symmetry breaking 
in 3 -I- 1 dimensions. The right graph of Fig. 6 shows the evolution for a 
sufficiently negative and a positive initial mass parameter M|, which leads 
to a nonzero and a vanishing field expectation value respectively. 

5. Applications 

5.1. Parametric resonance in quantum field theory 

In classical mechanics parametric resonance is the phenomenon of resonant 
amplification of the amplitude of an oscillator having a time-dependent 
periodic frequency. In the context of quantum field theory a similar phe- 
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Figure 7. Total energy (solid line) and classical field energy (dotted line) as a function of 
time for A = 10~®. The dashed line represents the fluctuation part, showing a transition 
between a classical field and a fluctuation dominated regime. 

nomcnon describes the amplification of quantum fluctuations, which can be 
interpreted as particle production. It provides an important building block 
for our understanding of the (pre)heating of the early universe at the end 
of an inflationary period^^, and may also be operative in various situations 
in the context of relativistic heavy- ion collisions'^. 

The phenomenon of parametric resonance provides a paradigm for the 
dynamics of quantum fields at nonperturbatively large densities'^. It is a 
far-from-equilibrium phenomenon involving the production of particles with 
densities inversely proportional to the coupling. For this reason, it has been 
much studied in the classical field approximation, which has long been the 
only available quantitative approach'^. Studies in quantum field theory 
have so far been limited to linear or mean- field type approximations'^. 

Using the methods described in the previous sections, recently the first 
quantum field theoretical study of this phenomenon beyond LO or Hartree 
approximations has been performed^^. To deal with the nonperturbatively 
large densities, one can employ the 2PI 1/A^-expansion at NLO, in the 
presence of a nonvanishing field (the "inflaton" ) in 3 + 1 dimensions. 

We consider"'^'^ a system that is initially in a pure quantum state, char- 
acterized by a large field amplitude (/)a = 0) ~ 1/ a/A and comparatively 
small quantum fluctuations, described here by employing vanishing parti- 
cle numbers at initial time. Fig. 7 provides an overview of the dynamics: 
At early times, the total energy of the system is dominated by the con- 
tribution from the large field amplitude, which subsequently decays into 



"^Another important example, where nonperturbatively high densities occur, is provided 
by the physics of the color glass condensate^® in QCD. 
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particles, giving rise to a fluctuation dominated regime at late times. More 
precisely, the early time coherent oscillations of the field trigger an expo- 
nential amplification of quantum fluctuations, corresponding to explosive 
particle production in a narrow range of momenta centered a,t p = po: this 
is parametric resonance. Later on, nonlinear interactions between field 
modes cause this amplification to propagate to higher momenta. More 
specifically, the initially amplified modes act as a source for other modes. 
This source-induced amplification results in an enhanced particle produc- 
tion in a broad momentum range. This is illustrated in Fig. 8, where the 
efi'ective particle number is shown for various momenta as a function of 
time from a numerical solution. 

We emphasize that, using the above 2PI techniques, one can obtain 
the main features of the nonlinear dynamics up to inonpert also from an 
analytical solution of the evolution equations^^. In particular, in order 
for the NLO efl'ects to be small compared to LO dynamics at early times, 
one needs values of A'' as large as A'^ > 1/A. The late-time dynamics for 
t ^ inonpert requires the numerical solution of the NLO equations. For 
the employed small couplings, this involves extremely large times which 
has not been undertaken up to now. It has recently been shown using 
classical field theory methods that the very slow subsequent evolution is 
characterized by turbulent behavior^°. However, to describe the late-time 
approach to quantum thermal equilibrium with a Bose-Einstein distributed 
particle number, one needs to go beyond the classical field approximation. 




20 40 60 80 



Figure 8. Effective particle number density for the transverse modes as a function of 
time for various momenta < p < 5po and A = 10~^. At early times, modes with 
momentum p ~ po are exponentially amplified with a rate ~ 270 via the mechanism 
of parametric resonance. Due to nonlinearities, one observes subsequently an enhanced 
growth with rate ~ 670 for a broad momentum range. 



12 



5.2. Nonequilibrium chiral quark-meson model 

Wc consider a theory involving two Dirac fcrmion flavors ( "quarks" ) coupled 
to a scalar a-field and a triplet of pseudoscalar "pions" as described by 
Eqs. (6). Such type of linear fr-models, which implement the approximate 
chiral S'?7l(2) x SUji{2) symmetry of QCD, can be a useful first guide for 
the behavior of strongly interacting matter relevant for the dynamics of 
relativistic heavy-ion collisions. 

Recently^^ the far-from-equilibrium time evolution and subsequent ther- 
malization has been described for this model at lowest nontrivial order in 
the 2PI coupling-expansion, corresponding to the two-loop graph presented 
in Fig. 3. Starting from various different far-from-equilibrium initial con- 
ditions, one finds a universal late-time behavior that is only determined by 
the expectation value of the initial energy density. In particular, one is able 
to observe the approach to Bose-Einstein and Fermi-Dirac distributions at 
sufficiently late times. Similar studies of the late-time behavior have previ- 
ously only been performed for purely scalar theories in 1 -f 1 dimensions''^. 

As an example, in Figs. 9 and 10 the statistical two-point functions 
for the fermions, Fv{t,t;p), and for the scalars, Frj,{t,t;p), are shown for 
two very different nonequilibrium initial conditions (cf. Ref. 11 for details). 
One observes that the time evolution becomes soon rather insensitive to 
the details of the initial condition. The time for the effective loss of initial 
conditions is well described by the inverse damping rate obtained from the 
respective unequal-time two-point function (cf. the example in Fig. 4). In 

0.3 
0.25 
0.2 

^ 0.15 

0.05 


-0.05 

5 10 15 20 25 30 50 100 500 

t 

Figure 9. Time evolution of fermionic equal-time two-point functions for various mo- 
menta and for two different initial conditions A and B. The corresponding initial particle 
number distributions, shown in the insets, are chosen such that the total energy is the 
same for both cases. (In units of the scalar thermal mass m.) 
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contrast, this time scale does not characterize the late-time behavior. For 
the latter, one finds, to very good approximation, an exponential relaxation 
of each mode to their universal late-time values. For the employed initial 
conditions the thermalization rate is somewhat larger than the damping 
rate. Though quantitatively very different, a similar qualitative behavior 
has been previously observed for scalar theories in 1 + 1 dimensions''^. 

At sufficiently late times one can explicitly demonstrate the approach to 
quantum thermal equilibrium*^ : Out of equilibrium, the spectral and statis- 
tical two-point functions are completely independent in general. However, 
if thermal equilibrium is approached, they have to become related by the 
fluctuation-dissipation relation at late times as stated in Eq. (14). For suf- 
ficiently late times one observes that the correlators become approximately 
homogeneous in time and a Fourier transform with respect to t — t' can be 
performed to very good approximation. The result is shown in Fig. 11. 

This work provides an important first step for a quantitative description 
of realistic theories with fermions in 3 -|- 1 dimensions, with important phe- 
nomenological applications^^ to the out-of-equilibrium chiral phase transi- 
tion, or the physics of baryogenesis in the early universe^^. 



We emphasize that thermal equilibrium cannot be reached on a fundamental level 
from time-reversal invariant evolution equations at any finite time. The results 

demonstratc^'^'^^ that thermal equilibrium ean be approached very closely at sufficiently 
late time, without again deviating from it for practically accessible times. 
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Figure 11. The late-time ratio of the statistical two-point function and the spectral func- 
tion in frequency space, both for fermions {Fy / py) and for scalars {F^/ p^). For thermal 
equilibrium the quotient has to correspond to the Bose-Einstein (BE) distribution for 
scalars and to the Fermi-Dirac (FD) distribution for fermions. The BE/FD distributions 
are displayed by the continuous curves with the same temperature T = 0.94m. 

6. Outlook 

The 2PI generating functional for Green's functions provides a powerful 
technique for nonequihbrium physics. Systematic approximation schemes 
can describe far-from-equilibrium dynamics as well as subsequent thermal- 
ization from "first principles" . Substantial progress has been achieved along 
these lines in recent years and, by now, the methods provide a practical 
means for quantitative descriptions of realistic scalar and fermionic quan- 
tum field theories. Apart from the discussed applications motivated by 
high-energy particle physics and cosmology, we emphasize that the same 
techniques can be applied to condensed matter systems. As an important 
example, the dynamics of Bose-Einstein condensation can be addressed 
along very similar lines. 

One important issue, which we have not discussed here, concerns the 
renormalizability of approximations based on truncations of the 2PI effec- 
tive action. To prove this is a nontrivial task because of the fully self- 
consistent character of these approximations. Recent work has shown that, 
in the context of scalar field theory, approximations based on a systematic 
loop-expansion of the 2PI effective action ( "^-derivable approximations") 
can be renormalized order by order, if the underlying theory is perturba- 
tively renormalizable^^'^^'^^. 

Another, major open question concerns the description of gauge fields. 
The difhculty is related to the fact that truncations of the 2PI effective 



15 



action in general violate Ward identities. One possibility to avoid this diffi- 
culty is to modify the 2PI-based approximations in such a way as to enforce 
Ward identities. Practical examples in this direction have been pointed out 
for nonequilibrium systems'^^. Another interesting possibility is to try to 
keep the gauge dependence of physical results under control. This has re- 
cently been investigated in Ref. 37 and further interesting developments are 
to be awaited. 
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